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SI.  Further  Background  on  Simulation  Methods 

In  summary,  we  employ  a  classical  atomistic  molecular  dynamics  approach  to  analyze 
specific  peptide  binding  to  metal  surfaces  in  comparison  with  thermochemical,  IR,  NMR,  and 
TEM  measurements  (Figure  1).  The  availability  of  suitable  parameters  for  biomolecules  and 
recent  improvements  in  the  accuracy  of  parameters  for  inorganic  components  by  more  than  one 
order  of  magnitude  ’  ~  are  thereby  very  helpful.  ^  Moreover,  modeling  at  all  levels, 
molecular,  coarse-grain,  and  bioinformatics,  is  ultimately  important  in  identifying  the  most 
suitable  peptide  sequences  for  controlled  binding  and  detachment  (Figure  1)  as  the  huge  number 
of  possible  sequences  exceeds  the  number  of  available  sequences  in  commercial  phage  libraries 
by  more  than  a  million.  Quantum-mechanical  approaches  indicate  approximate  trends  of  the 
interaction  of  peptide  fragments  and  a  few  solvent  molecules  with  parts  of  a  surface,  and  have 
shown  that  covalent  interactions  with  metal  surfaces  are  modest  to  small.  >  •  •  ’  •  However, 
such  approaches  cannot  fully  explain  the  mechanisms  of  binding  due  to  limitations  to  static 
calculations  including  a  few  hundred  atoms,  and  neglect  of  thousands  of  needed  solvent 
molecules,  ions,  and  dynamics  at  significant  time  scales.  Some  ab-initio  results  also  disagree 
with  each  other  due  to  known  difficulties  with  the  treatment  of  van-der-Waals  interactions  (only 
in  wave-function  based  techniques)  and  accurate  electrostatics. A  limitation  of  ab-initio 
and  classical  atomistic  methods  is  also  the  difficulty  to  include  polarization  effects  on  metal 
surfaces though  estimates  suggest  modest  contributions  on  even  surfaces.  Coarse-grain 
approaches  can  be  computationally  ~10  times  more  efficient  and  include  certain  specific 
peptide-surface  interactions.  A  critical  role  can  also  be  attributed  to  bioinformatics  approaches 
such  as  simple  numerical  screening  functions  on  the  basis  of  molecular-level  insight  to  help 
eliminate  sequences  of  undesirable  binding  strength.  The  success  of  such  higher-level 
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approaches,  however,  yet  depends  on  understanding  the  nature  of  the  molecular  interactions 
(Figure  1). 

S2.  Efficient  Computation  of  Adsorption  Energies 

The  quantitative  computation  of  adsorption  energies  requires  the  inclusion  of  explicit 
water  molecules  and  ionic  solutes  due  to  the  formation  of  spatially  oriented  hydrogen  bonds 
which  change  their  orientation  when  peptides  (or  other  adsorbates)  move  from  the  solution  to  the 
surface.  The  calculation  of  adsorption  energies  also  requires  an  accurate  summation  of  Coulomb 
interactions  on  the  order  of  <0.1  kcal/mol  uncertainty  for  the  entire  system  (see  section  S3. 3). 
Under  these  conditions,  we  have  developed  and  tested  two  methods  to  compute  the  changes  in 
energy,  entropy,  and  free  energy  upon  adsorption  (Figure  SI). 

S2.1.  Method  1  -  Simple  Approach.  A  conceptually  simple  approach  involves  the 
simulation  of  the  adsorbate  on  the  surface  (simulation  1)  and  in  solution  (simulation  2)  using  the 
same  box  size  and  box  content  (Figure  SI  a).  The  box  needs  to  be  large  enough  to  avoid  contact 
of  the  adsorbate  with  the  surface  during  the  solution  run.  However,  large  boxes  result  in  high 
values  and  large  fluctuations  of  the  total  energy  which  require  extra  simulation  time  to  compute 
the  relatively  small  energy  difference  AE  =  E^-E^  between  simulation  1  and  simulation  2  with 
statistically  acceptable  errors  (Figure  SI  a).  The  time  dependence  on  the  number  of  particles  N  is 
typically  on  the  order  N  for  summations  of  Coulomb  interactions  using  the  Ewald  method  and 
for  summations  of  van-der-Waals  interactions  using  spherical  cutoffs,  and  on  the  order  N-\r\N  for 
summation  of  Coulomb  interactions  using  the  PPPM  method. 

Method  1  is  less  efficient  due  to  the  inclusion  of  many  water  molecules  in  the  surface 
run,  and  due  to  the  inclusion  of  the  surface  atoms  and  surrounding  water  molecules  in  the 
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solution  run  (Figure  Sib).  This  introduces  particular  computational  overhead  when  the 
adsorption  of  a  series  of  molecules  on  the  same  surface  is  of  interest.  In  the  presence  of  usual 
periodic  boundary  conditions,  these  issues  of  method  1  can  be  circumvented  by  taking  the 
additivity  of  intermolecular  interactions  into  account.  As  shown  earlier, the  interaction  of  polar 
but  overall  charge  neutral  molecules  or  crystals  becomes  almost  negligible  at  distances  greater 
than  3  nm  so  that  larger  simulation  boxes  in  Figure  Sla  can  be  split  into  multiple  boxes  of  ~5  nm 
height  as  shown  in  Figure  Sib.  This  partition  of  the  computation  boxes  leads  to  method  2. 

S2.2.  Method  2  -  Efficient  Approach.  In  the  second  scheme  (Figure  Sic),  a  total  of  four 
calculations  with  smaller  boxes  are  required  which  is  approximately  twice  as  fast  as  two 
calculations  with  larger  boxes  using  Ewald  procedures  or  spherical  cutoffs  for  Lennard-Jones 
interactions  (t(N)  ~  N^).  For  computational  screening  of  additional  adsorbates  on  the  same 
surface,  only  two  additional  calculations  per  adsorbate  are  required  to  obtain  Esutf+A+Hio  and 
Ea+h2o  while  the  values  for  Ehio  and  Ea+h2o  remain  the  same.  Assuming  a  dependence  of 
simulation  time  t(N)  on  the  number  of  atoms  N,  a  box  size  of  2N  in  Fig.  Sla,  and  about  half  the 
box  size  N  in  Fig.  Sic,  the  computation  of  adsorption  energies  for  k  peptides  on  the  same  surface 


takes  the  relative  time 


2k +  2 
2k 


t(N) 

t(2N) 


For  example,  method  2  will  screen  20  peptides  in  26%  of 


the  time  compared  to  method  1  in  a  simulation  with  t(N)  ~  (Ewald  or  spherical  cutoffs),  and 
in  less  than  50%  of  the  time  compared  to  method  1  in  a  simulation  with  t(N)  ~  N  In  N  (PPPM 
summations). 

A  further  advantage  of  method  2  over  method  1  is  the  reduction  of  the  total  energy  and  its 
standard  deviation,  which  lowers  the  standard  deviation  of  computed  adsorption  energies  on  an 
absolute  scale.  Possible  interferences  from  the  surface  in  the  solution  run  in  method  1  are  also 
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eliminated,  and  method  2  equally  provides  access  to  thermodynamic  quantities  including 
energies,  entropies,  and  free  energies  of  adsorption.  The  computation  of  entropy  AS  (using 
A  =  U  -TS  or  other  method)  and  free  energy  AA  is  more  time-consuming  than  that  of  the 
energies,  however,  because  very  well  resolved  histograms  of  the  energy  distribution  are  needed 
or  multiple  simulations  for  each  system  to  facilitate  thermodynamic  integration,  respectively. 

The  application  of  either  method,  provided  the  box  size  is  large  enough  to  eliminate  finite 
size  effects  and  simulations  are  allowed  to  reach  equilibrium,  leads  to  identical  results. 
Simulations  are  best  performed  using  the  NVT  ensemble  and  additive  molecular  volumina,  i.e., 
predefined  densities  for  each  constituent  such  as  metal  atoms,  water  molecules,  and  peptides. 
Then,  for  example,  using  the  same  cross  section  (A  =  xy)  of  the  box,  the  box  height  (z)  is 
incrementally  additive  for  every  component  added.  To  compute  differences  in  thermodynamic 
quantities,  the  simulation  temperature  must  be  constant  or  corrected  to  the  same  value  using  the 
heat  capacity  of  each  box  (as  derived  in  the  simulation). 
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(a)  Method  1 


(b)  Decomposition  as 
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(c)  Method  2 


AE=  Esurf+A+H20  ~  Ea+H20  +  ^H20  ”  ^Surf+H20 


AS,  AG  analogous 
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Fig.  SI.  Computation  of  adsorption  energies  under  periodic  boundary  conditions  (illustrations  in  2D  for  simplicity),  (a)  Method  1 
involves  the  calculation  of  average  energies  in  two  simulation  boxes  with  the  adsorbate  on  the  surface  and  in  solution,  (b) 
Decomposition  into  smaller  boxes  for  reduced  computation  time  and  lower  fluctuation  in  computed  energy  differences,  (c)  Method  2 
involves  four  smaller  boxes  and  is  more  efficient  than  method  1,  particularly  for  screening  many  molecules  on  the  same  surface. 
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S3.  Computational  Details 

S3.1.  Models  of  the  Peptides,  Water,  and  the  Metallic  Substrates.  Atomistic  Models 
of  the  peptides  were  build  from  scratch  using  the  amino  acid  sequences  and  protonation  state  at 
pH  =  7  outlined  in  Table  1.  In  addition  to  the  possible  presence  of  counterions  (Flg-Naa,  Flgd- 
Nai,  Pd2-Cl,  Pd4-Cl),  all  peptides  are  used  in  their  zwitterionic  form.  Models  of  the  peptides  in 
different  backbone  conformations  (a-helix,  extended,  random  coil)  were  prepared  with  the  help 
of  the  Hyperchem  visualizer,"^^  as  well  as  by  manual  drawings  (random  coil)  using  the  Materials 
Studio  builder  and  simple  energy  minimization.  For  water,  we  employ  the  standard  SPC  model 
implemented  in  the  consistent  valence  force  field  (CVFF)."^^ 

Models  of  the  fee  metals  were  constructed  from  the  unit  cell  parameter  Uq  and  the  atom 

positions  (0,0,0),  (1/2,172,0),  (1/2,  0,  1/2),  (0,  1/2,  1/2).^^  The  fee  unit  cell  and  related  supercells 
yield  {100}  surfaces  along  the  Cartesian  coordinate  axes.  To  construct  {111}  surfaces,  an 


equivalent  larger  cell  was  build  as  ^ 


:  VSag  with  the  atom  positions  (0,  0,  0), 


(1/2,  1/2,  0),  (0,  1/3,  1/3),  (1/2,  5/6,  1/3),  (1/2,  1/6,  2/3),  (0,  2/3,  2/3).  {Ill}  surfaces  are  then 
obtained  perpendicular  to  the  z  axis.  The  Pd-Au  {111}  bimetal  interface  was  constructed 
starting  with  the  Au  supercell  and  replacement  of  half  the  atoms  by  Pd.  The  lattice  constants  are 

o  o 

very  similar  (4.0782  A  vs  3.8903  A)  so  that  the  scale  of  the  simulation  box  (3  nm)  does  not  lead 
to  dislocations;  a  fractional  offset  of  the  Pd  lattice  positions  from  the  nominal  Au  lattice 
positions  was  observed  in  the  course  of  the  simulation. 

Complete  start  structures  consist  of  the  inorganic  substrate,  the  peptide,  and  water,  or 
subsets  thereof  with  three-dimensional  periodic  boundary  conditions  (see  Figure  Sic).  The 
structures  for  each  set  of  simulations  were  assembled  from  slabs  of  the  metal  substrate,  pre¬ 
equilibrated  assemblies  of  1000  water  molecules  (2000  water  molecules  using  method  1  in 
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Figure  Sla),  and  one  peptide  chain  using  the  Materials  Studio  graphical  interface.  The  initial 
conformations  of  the  peptide  chains  were  obtained  from  the  build  of  various  conformations,  the 
result  of  the  solution  run,  as  well  as  from  pre-equilibrated  peptide  conformation  on  the  metal 
surface  in  vacuum  (without  solvent).  The  setup  of  the  hybrid  structures  assumes  a  density  of 
1000  kg/m  for  all  liquid  phases  such  as  water  and  peptide  solutions,  and  the  exact  density  of  the 
metal  consistent  with  unit  cell  parameters.  The  box  dimensions  are  typically  3  nm  x  3  nm  x  z 
whereby  the  box  height  z  amounts  to  >1.2  nm  for  metal  substrates,  >3.5  nm  for  water,  and  >3.5 
nm  for  the  aqueous  solutions  of  the  peptides.  The  size  was  chosen  such  that  25%  to  40%  of  the 
surface  area  is  covered  by  the  peptide,  close  to  experimental  conditions.  This  comparatively  low 
surface  coverage  for  a  single  peptide  chain  in  the  simulation  box  also  reduces  finite  size  effects 
to  a  tolerable  level.  The  peptides  can  be  initially  superimposed  onto  the  water  slab;  subsequent 
energy  minimization  removes  close  contacts  with  only  modest  distortions  of  the  initial  secondary 
structure.  The  concept  of  additive  molecular  volumes  and  the  choice  of  a  benchmark  density  for 
every  component  is  required  for  NVT  simulations  according  to  Figure  S 1 ;  computed  adsorption 
properties  are  not  affected  if  benchmark  densities  are  changed  consistently  in  all  simulations  in  a 
range  of  +1%. 

S3.2.  Force  Field  Parameters  for  Peptides  and  Metals.  At  the  classical  atomistic  level, 
several  force  fields  with  good  parameters  for  peptides  and  water  are  available,  e.g.,  AMBER, 
CHARMM,  CVFF,  OPLS-AA,  PCFF,  and  compatible  parameters  for  inorganic  structures  have 
been  recently  developed.  ’  The  improvement  of  Lennard-Jones  parameters  for  fee  metals 
leads  to  better  than  10%  agreement  of  surface  and  interface  properties  with  experiment  for  fee 

OO  'iO  'iA 

metals  ’  ~  which  is  essential  to  understand  sensitive  interfacial  adsorption  processes.  In 
contrast,  earlier  metal  parameters  in  force  fields  were  associated  with  deviations  in  excess  of 
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100%  and  result  in  wide  gaps  between  computed  adsorption  properties  and  actual  interfacial 

3  23—27 

interactions  in  experiment.  ’  In  this  paper,  we  employ  the  consistent  valence  force  field 
(CVFF)"^*’"^^  extended  for  accurate  Lennard-Jones  parameters  for  fee  metals  (Ag,  Al,  An,  Cu,  Ni, 
Pb,  Pd,  Pt).^*’ 


This  extended  version  of  CVFF,  therefore,  contains  the  existing  parameters  for  the 
peptides  and  water"'^^’"^^  as  well  as  new,  accurate  parameters  for  fee  metals. The  potential  energy 
expression  with  no  cross-terms  and  no  morse  potential  is  employed: 

F  =F  +  F  4-F  +F  F  +  F  tSIl 

pot  ^  bond  ^  angle  ^torsion  ^  out— of -plane  Coulomb  ^vdW  ^ 


“  ^0,y)  ■*"  ^O.yH  )] 

ij  bonded  ijk  bonded  ijkl  bonded 


+  Eoop  + 
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A  'TTP  p  Y 

“r/tOQC.^  // nonbonded  '  ij  tynonbonded 


(1,2  and  l,3excl) 


(1,2  and  1, 3excl)  V 


'ojj 

J 


-2 


r  r  ^ 

'o,ij 


6  A 


J 


Equation  (SI)  shows  the  additive  energy  terms  for  quadratic  bond  stretching  ,  quadratic 


angle  bending  ,  a  one-term  torsion  potential  ,  quadratic  out-of-plane  deformation 


Ecu, -of -planet  clectrostatic  energy ,  and  a  12-6  Lennard-Jones  potential  for  van-der- 

Waals  interactions.  Nonbond  interactions  between  1,4  covalently  bonded  atoms  are  included  to 
100%  and  geometric  combination  rules  apply  for  Lennard-Jones  parameters  between  different 
atom  types.  A  critical  advantage  over  earlier  models  is  the  accurate  12-6  Lennard-Jones  potential 
for  fee  metals.  Computed  cell  parameters  for  the  metals  approach  experimental  values  within 
0.2%  deviation,  computed  surface  tensions  for  the  {111}  surfaces  of  the  metal  fall  within  1-5% 
of  the  experimental  values  at  298  K  (equal  to  the  uncertainty  in  experiment),  and  the  computed 
surface  energy  anisotropy  between  the  {111}  and  {100}  metal  surfaces  compares  favorably  with 
experimental  data.  ’  ’  Computed  metal- water  interfacial  tensions  are  on  average  10%  to  15% 
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lower  than  expected  from  experimental  data  according  to  the  Young  equation.  A  lower 
interface  tension  might  have  some  justification,  however,  due  to  attractive  polarization  effects  as 
the  Young  equation  only  provides  an  upper  bound.  In  detail,  polarization  will  be  discussed  in  a 
separate  contribution.  In  summary,  computed  surface  energies  using  metal  parameters  in  earlier 
semiempirical  force  fields  were  associated  with  errors  on  the  order  of  100%  so  that  the  present 
semi-quantitative  treatment  represents  an  improvement  of  between  one  and  two  orders  of 
magnitude  in  accuracy.  Remaining  sources  of  uncertainties  are  (1)  polarization  effects,  (2)  the 
approximate  nature  of  the  energy  model,  specifically  limitations  of  the  SPC  water  model  and  of 
the  peptide  parameters,  (3)  geometric  combination  rules  of  the  12-6  Lennard-Jones  parameters, 
(4)  and  the  inability  of  the  force  field  to  consider  covalent  contributions,  such  as  major  shifts  in 
electron  density,  to  adsorption. 

S3.3.  Simulation  Protocol.  The  programs  Discover'^^  and  LAMMPS^'^  were  used  to  carry 
out  molecular  dynamics  (MD)  simulations.  Typically,  we  started  with  a  short  energy 
minimization  for  200  steps  using  the  conjugate  gradient  method,  and  then  carried  out  NVT 
molecular  dynamics.  A  time  step  of  1  fs,  the  Verlet  integrator,  the  Anderson  or  the  Nose 
thermostat  at  298.15  K,  a  van-der-Waals  cutoff  at  1.2  nm,  Ewald  summation  of  Coulomb 
interactions  with  high  accuracy  (10“^  kcal/mol),  or  PPPM  summation  of  Coulomb  interactions 
with  high  accuracy  (10“^  kcal/mol)  were  employed.  The  analysis  of  adsorption  energies  and 
chain  conformations  followed  as  a  post-processing  operation  after  typical  simulation  times  of  5 
ns  to  10  ns,  corresponding  to  5  to  10  million  time  steps.  For  the  analysis,  the  first  1  ns  or  more  of 
the  trajectory  was  discarded  to  take  into  account  only  the  parts  of  the  trajectory  with  steady 
energies  and  thermodynamically  significant  conformations.  We  note  thermodynamic  averages. 


11  of  16 


as  well  as  the  visual  and  statistical  inspection  of  complete  trajectories  are  required  to  understand 
conformational  changes  so  that  representative  snapshots  are  intended  for  guidance  only. 

A  bottleneck  is  the  generation  of  a  Boltzmann  average  of  equilibrium  conformations  of 
the  peptides  on  the  surface.  While  the  modest  length  of  the  peptides  up  to  12  amino  acids  does 
not  excessively  complicate  this  goal,  several  start  conformations  that  are  distinct  from  each 
other,  such  as  helix,  outstretched,  random  coil,  and  optimized  structures  in  vacuum  were 
required.  Solution  structures  converge  to  consistent  results,  however,  on  the  metal  surfaces  with 
particularly  strong  adhesion,  we  cannot  be  certain  that  the  global  minimum  energy  has  been 
reached  in  every  case  as  relaxation  times  exceed  10  ns.  However,  the  deviation  in  adsorption 
energies  is  mostly  in  a  range  of  +5  kcal/mol  for  different  start  structures  in  independent 
simulations  which  is  less  than  10%  of  the  strongest  adsorption  energies  seen.  Besides,  reference 
to  coarse-grain  models  and  Monte  Carlo  algorithms  to  pre-equilibrate  the  structures  in  a 
sequential,  reversible  two-scale  approach  may  not  result  in  major  improvements  in  sampling  as 
the  competition  of  water  and  peptide  on  the  surface  and  specific  polar  chemistries  of  side  groups 
are  difficult  to  capture  in  simplified  models. 

An  important  aspect  in  the  computation  of  adsorption  energies  was  the  adjustment  of  the 
average  total  energy  for  each  simulation  to  the  same  reference  temperature  (298.15  K)  using  the 
heat  capacity  of  the  simulation  box.  Even  small  temperature  differences  (<0.1  K)  between 
individual  calculations  (Figure  SI)  cause  noticeable  energy  differences  (<2  kcal/mol)  and  could 
introduce  errors  in  computed  adsorption  energies. 
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S4.  Additional  Evidence  on  the  Binding  Mechanism  from  Previous  Simulations 

Additional  evidence  for  the  proposed  binding  mechanism  by  earlier  computational  results  at 
various  length  scales  is  described  in  the  following. 

S4.1.  Ab-initio  methods.  Ab-initio  methods  have  been  employed  to  study  peptide 
fragments  such  as  amino  acids  or  water  in  contact  with  fee  metals  at  the  smallest  level  without 
significant  dynamics  or  solvation.  ’  ’  ’  ’  In  a  benchmark  study,  Haftel  et  al.  compared  the 
reliability  of  tight-binding  (TB),  electron  density  functional  methods  (DFT),  embedded  atom 
models  (EAM),  and  modified  embedded  atom  models  (MEAM)  to  compute  surface  tensions  of 
several  fee  metals  (Rh,  Ir,  Ni,  Pd,  Pt,  Au).  The  methods  exhibit  scatter  by  a  factor  of  two  among 
each  other,  and  typical  deviations  from  experiment  are  -20%  in  each  method  with  individual 
deviations  up  to  50%  both  above  and  below  experimental  values.  The  reliability  is  therefore  a 
problem  and  must  be  taken  into  account  in  the  discussion  of  DPT  results,  as  well  as  difficulties  to 
combine  DPT  methods  with  classical  force  fields  to  simulate  the  dynamics  of  peptides  with 
explicit  solvent  molecules  in  solution. 

Water  adsorption  and  layering  on  transition  metal  surfaces  was  studied  by  Vassilev  et 
al.  using  Born-Oppenheimer  DPT.  On  Rh  {111}  and  Pt  {111}  surfaces,  a  planar  orientation  of 
water  molecules  is  found,  and  the  adsorption  energy  of  single  water  molecules  in  vacuum  at  low 
surface  coverage  was  computed  as  8.6  and  6.9  kcal/mol,  at  a  higher  coverage  with  a  water 
bilayer,  13.1  kcal/mol  and  13.1  kcal/mol  for  Rh  and  Pt  surfaces.  Schravendijk  et  al.  also  found 
a  planar  coordination  of  water  molecules  on  the  metal  surfaces  and  computed  slightly  higher 
adsorption  energies  of  single  water  molecules  of  9.7  and  8.1  kcal/mol  on  Rh  {111}  and  Pt  {111} 
surfaces,  as  well  as  2.3  and  7.6  kcal/mol  on  Au  {111}  and  Pd  {111}  surfaces.  Classical  MD 
simulation  (section  S3)  agrees  with  the  preference  for  a  planar  orientation  and  adsorption 
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energies  of  4.4  and  5.4  kcal/mol  are  computed  on  the  An  {111}  and  on  the  Pd  {111}  surface  at 
<5%  surface  coverage  (calculation  as  in  Figure  Sla  with  one  water  molecule  as  adsorbate  in 
vacuum).  These  values  agree  with  lower  surface  tensions  of  Au  and  Pd  versus  Rh  and  Pt;  the 
gap  in  DFT  calculations  between  the  Au  {111}  and  Pd  {111}  surface  yet  remains  to  be  explained 
on  physical  grounds.  Experimental  measurements  of  contact  angles  indicate  0°  and,  therefore, 
hydrophilic  properties  of  both  clean  Au  and  Pd  surfaces,^^  in  agreement  with  attraction  of  water 
to  the  metal  surfaces  in  the  computation.  Further,  dual-scale  simulations  of  benzene  adsorption 
on  Ni  {111}  and  Au  {111}  surfaces  including  a  few  explicit  water  molecules  were  carried  out  by 
DFT  and  Car-Parrinello  ab-initio  MD  molecular  dynamics  as  well  as  specialized  Lennard-Jones 
and  Morse  potentials  for  a  classical  description  of  the  interaction  with  the  metal.  Taking  into 
account  competition  between  a  few  water  molecules  and  the  benzene  solute,  the  adsorption 
energy  of  benzene  in  solution  on  the  Ni  {111}  surface  was  computed  as  9.6  kcal/mol,  slightly 
stronger  adsorption  was  found  on  the  Pd  {111}  surface,  and  a  preference  for  desorption  was 
predicted  on  the  Au  {111}  surface.  However,  the  surface  tension  of  the  metals  was  not  validated 
in  these  models  and  the  few  included  water  molecules  are  a  crude  approximation  of  a  solvation 
shell.  Nevertheless,  the  binding  energy  agrees  with  the  order  of  magnitude  for  aromatic  amino 
acids  (Table  2  and  Table  4),  and  a  weaker  attraction  of  benzene  (in  Phe)  to  the  Au  {111}  surface 
is  also  observed  in  our  study.  Our  results  indicate  still  adsorption  in  agreement  with  experimental 
data  (section  5.1).  ’  ’  It  may  be  alternatively  possible  that  phenyl  rings  are  positioned  flat  on  the 
surface  only  in  order  to  provide  sufficient  “contact  area”  for  neighboring  parts  of  the  molecule 
which  are  factually  involved  in  an  adsorption  process.  Schravendijk’s  method  was  further 
extended  into  a  highly  specialized  force  field  description  for  the  interaction  of  selected  hydrated 
amino  acids  with  a  Ni  {111}  surface. A  new  extensive  set  of  parameters  would  be  required  for 
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every  additional  amino  acid  and  every  additional  metal.  Meanwhile,  the  validation  of  the  metal 
surface  energy  and  the  physical  justification  of  the  numerous  parameters  remain  unclear,  and 
incompatibility  with  existing  biomolecular  force  renders  the  approach  less  practical.  Therefore, 
atomistic  simulation  with  extended  force  field  engines  as  presented  in  this  paper  yields  results  of 
the  same,  if  not  better,  quality,  and  is  simpler,  faster,  and  more  widely  applicable. 

Ghiringhelli  et  al.  studied  Phe  adsorption  in  vacuum  on  {111}  surfaces  of  fee  metals 
using  DFT.  Strong  adsorption  of  -30  to  -25  kcal/mol  (Pt,  Ni,  Pd)  and  weaker  adsorption  of  -12 
to  -7  kcal/mol  (Cu,  Ag,  Au)  was  found.  For  Cu,  Ag,  Au,  the  phenyl  ring  is  not  in  close  contact 
with  the  surface.  Although  water  was  absent  in  this  study,  the  range  of  values  is  similar  to  those 
inferred  from  the  force-field  based  simulation  (Tables  2  and  Table  4);  a  lesser  attraction  of  the 
phenyl  ring  to  the  Au  {111}  surface  compared  to  the  Pd  {111}  surface  is  observed  in  our  method 
as  well.  Hong  et  al.  have  shown  by  further  DFT  calculations  in  vacuum  that  charge  transfer  is 
possible  on  Au  and  Pd  surfaces.  Up  to  -0.65e  could  be  transferred  per  negatively  charged 
residue  such  as  C02~  of  Asp  over  numerous  electron  accepting  superficial  Au  atoms  in  vacuum 
(Table  3).  In  aqueous  solution,  however,  the  presence  of  water  molecules  as  a  solvent  of  similar 
polarity  diminishes  the  amount  of  charge  transfer  from  a  single  COi”  group  which  would  be 
more  challenging  to  analyze  using  this  method.  Among  Asp,  Lys,  Arg,  Ser,  Pro,  and  Val,  the 
strongest  binding  residue  in  vacuum  has  been  Asp  by  a  considerable  margin  (-57  kcal/mol), 
followed  by  Lys  (-27  kcal/mol)  and  Arg  (-19  kcal/mol),  and  almost  no  binding  was  found  for 
Pro  and  Ser.  While  the  values  will  be  substantially  smaller  in  aqueous  solution,  this  ranking 
coincides  with  the  force-field  based  simulation  (Tables  2  and  4). 

S4.2.  Force-field  based  simulations.  Sarikaya  et  al.^’^^  carried  out  classical  molecular 
dynamics  simulations  of  large  proteins  (>40  amino  acids)  on  Au  {111}  and  {112}  surfaces  with 
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a  number  of  heavy  approximations  which  weigh  on  the  reliability  of  the  results:  (1)  The  proteins 
were  placed  in  a  dielectric  continuum  so  that  explicit  water  molecules  and  hydrogen  bonds  are 
absent.  (2)  The  surface  properties  of  Au  according  to  the  chosen  Lennard-Jones  parameters  were 
not  validated  which  adds  uncertainty  in  excess  of  100%  to  adsorption  energies  and  affects  the 
molecular  mechanism  of  adsorption  of  the  peptide.  (3)  The  thickness  of  the  surfaces  was 

o 

restricted  to  only  3  atomic  layers  (6  A)  which  is  shorter  than  the  range  of  interactions  between 

o 

the  metal  atoms  (at  least  10-12  A).  Comparisons  relative  to  experiment  are  rather  difficult  under 
these  conditions,  and  further  simulations  of  peptides  on  metal  surfaces  using  this  approach  lead 
to  a  number  of  questionable  specific  results.  An  interesting  outcome,  however,  is  the  response 
of  peptides  to  gross  surface  features  such  as  ridges  in  {110}  surfaces.  A  subsequent  study  by 

-yc 

Kantarci  et  al.  extended  the  approach  using  an  all  atom  representation  of  peptides  and  water  on 
Pt  {111}  surfaces.  This  represents  a  more  realistic  model,  however,  the  Lennard-Jones 
parameters  for  the  metal  according  to  the  unmodified  CVFF  overestimate  the  Pt  surface  tension 

9Q 

by  101%  relative  to  experiment  and  modify  the  interfacial  structures  and  adsorption  processes; 
the  thin  metal  layer  also  introduced  residual  interactions  between  water  molecules  above  and 
below  the  metal  slab.  Under  these  severe  approximations,  Ser,  Thr,  Arg,  and  Pro  were  found  in 
proximity  to  the  metal  surface  (Pt).  Overall,  earlier  all  atom  models  may  not  have  been  reliable 
enough  for  quantitative  or  firm  qualitative  conclusions. 
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